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Aims. Twisted magnetic fields are frequently seen to emerge above the visible surface of the Sun. This emergence is 
usually associated with the rise of buoyant magnetic flux structures. Here we ask how magnetic fields from a turbulent 
large-scale dynamo appear above the surface if there is no magnetic buoyancy. 

Methods. The computational domain is split into two parts. In the lower part, which we refer to as the turbulence zone, 
the flow is driven by an assumed helical forcing function leading to dynamo action. Above this region, which we refer 
to as the exterior, a nearly force- free magnetic field is computed at each time step using the stress- and-relax method. 
Results. Twisted arcade- like field structures are found to emerge in the exterior above the turbulence zone. Strong current 
sheets tend to form above the neutral line, where the vertical field component vanishes. Time series of the magnetic 
field structure show recurrent plasmoid ejections. The degree to which the exterior field is force free is estimated as the 
ratio of the dot product of current density and magnetic field strength to their respective rms values. This ratio reaches 
values of up to 95% in the exterior. A weak outward flow is driven by the residual Lorentz force. 

Key words, magnetohydrodynamics (MHD) - turbulence - Sun: dynamo - stars: magnetic fields 



1. Introduction 

The magnetic field at the visible surface of the Sun is known 
to take the form of bipolar regions. Above these magnetic 
concentrations the field continues in an arch-like fashion. 
These formations appear usually as twisted loop-like struc- 
tures. These loops can be thought of as a continuation of 
more concentrated flux ropes in the bulk of the solar con- 
vection zone. However, this interpretation is problematic 
because we cannot be certain that the magnetic field in the 
Sun is generated in the form of flux ropes. Indeed, simu- 
lations of successful large-scale dynamos suggest that con- 
centrated tube-like structures are more typical of the early 
kinematic stage, but in the final nonlinear stage the field 
becomes more space-filling (Brandenburg 2005, Kapyla et 
al. 2008). The idea that the dynamics of such tubes is gov- 
erned by magnetic buoyancy is problematic too, because 
the solar convection zone is strongly stratified with concen- 
trated downdrafts and broader upwellings. This leads to ef- 
ficient downward pumping of magnetic field toward the bot- 
tom of the convection zone (Nordlund et al. 1992; Tobias 
et al. 1998). This downward pumping is generally found 
to dominate over magnetic buoyancy. The question then 
emerges whether magnetic buoyancy can still be invoked as 
the main mechanism for causing magnetic flux emergence 
at the solar surface. Another possible mechanism for the 
emergence of magnetic field at the solar surface might sim- 
ply be the relaxation of a strongly twisted magnetic field 
in the bulk of the convection zone. Twisted magnetic fields 
are produced by a large-scale dynamo mechanism that is 
generally believed to be the motor of solar activity (Parker 
1979). One such dynamo mechanism is the a effect that pro- 
duces a large-scale poloidal magnetic field from a toroidal 
one. However, this mechanism is known to produce mag- 



netic fields of opposite helicity (Seehafer 1996; Ji 1999). 
This magnetic helicity of opposite sign is an unwanted by- 
product, because it quenches the dynamo effect (Pouquet 
et al. 1976). A commonly discussed remedy is therefore to 
allow the helicity of small-scale field to leave the domain, 
possibly in the form of coronal mass ejections (Blackman 
& Brandenburg 2003). 

In order to study the emergence of helical magnetic 
fields from a dynamo, we consider a model that combines a 
direct simulation of a turbulent large-scale dynamo with a 
simple treatment of the evolution of nearly force-free mag- 
netic fields above the surface of the dynamo. An additional 
benefit of such a study is that it alleviates the need for 
adopting a boundary condition for the magnetic field at 
the top of the dynamo region. This is important, because it 
is known that the properties of the generated large-scale 
magnetic field strongly depend on boundary conditions. 
A common choice for the outer boundary condition is to 
assume that the magnetic field can be matched smoothly 
to a potential field. Such a condition is relatively easily 
implemented in calculations employing spherical harmonic 
functions (see, e.g., Krause & Radler 1980). A more real- 
istic boundary condition might be an extrapolation to a 
force-free magnetic field where the Lorentz force vanishes 
in the exterior. This means that the current density is pro- 
portional to the local magnetic field, but the constant of 
proportionality depends generally on the magnetic field it- 
self. This renders the magnetic boundary condition nonlin- 
ear, which is therefore not easy to implement. Moreover, 
a perfectly force-free magnetic field is not completely re- 
alistic either. Instead, we know that above the solar sur- 
face, magnetic fields drive flares and coronal mass ejections 
through the Lorentz force. A more comprehensive approach 
would be to include in the calculations the exterior regions 



2 



Warnecke & Brandenburg: Surface appearance of dynamo-generated fields 



above the solar or stellar surface. This can be computa- 
tionally prohibitive and a realistic treatment of that re- 
gion may not even be necessary. It may therefore make 
sense to look for simplifying alternatives. One possibility is 
therefore to attempt an iteration toward a nearly force-free 
magnetic field such that the field can deviate from a force- 
free state locally in regions where the field cannot easily be 
made force-free. This could be done by solving the induc- 
tion equation with an additional ambipolar diffusion term, 
which implies the presence of an effective velocity correc- 
tion proportional to the local Lorentz force. This is some- 
times called the magneto-frictional method and has been 
introduced by Yang et al. (1986) and Klimchuk & Sturrock 
(1992). In this approach the electromotive force attains not 
only a term in the direction of B, but also a term in the 
direction of J (Brandenburg & Zweibel 1994). The latter 
corresponds to a diffusion term, which explains the diffu- 
sive aspects of this effect. However, the resulting ambipolar 
diffusivity coefficient is proportional to B 2 and can locally 
become so large that the computational time step becomes 
significantly reduced. This is a typical problem of parabolic 
equations. A better method is therefore to turn the prob- 
lem into a hyperbolic one and to solve an additional evo- 
lution equation for the velocity correction where the driv- 
ing force is the Lorentz force. This approach is sometimes 
called the "force- free model" (FFM), even though the field 
in this model is never exactly force-free anywhere (Mikic 
et al. 1988; Ortolani & Schnack 1993). In the context of 
force-free magnetic field extrapolations this method is also 
known as the stress-and-relax method (Valori et al. 2005). 

2. Equations for the Force-Free Model 

The equation for the velocity correction in the Force-Free 
Model (FFM) is similar to the usual momentum equation, 
except that there is no pressure, gravity, or other driving 
forces on the right-hand side. Thus, we just have 

DU 



Dt 



J x B/p + F v 



(1) 



where J x B is the Lorentz force, J = V x B / po is the 
current density, po is the vacuum permeability, -F v isc is the 
viscous force, and p is here treated as a constant corre- 
sponding to a prescribed density. Equation (pQ) is solved 
together with the induction equation. In order to preserve 
V • B = 0, we write B = V x A in terms of the vector 
potential A and solve the induction equation in the form 
ft A 

^-=UxB + r ] V 2 A, (2) 

where we have adopted the so-called resistive gauge in 
which the electrostatic potential is equal to — t?V • A and 
the magnetic diffusivity r] is assumed constant. The value of 
rj will be given in terms of the magnetic Reynolds number, 
whose value will be specified below. No continuity equation 
for p is solved in this part of the domain, because there is 
no pressure gradient in the momentum equation. 

In the following we couple such a model for the magnetic 
field above the solar photosphere to a simulation of a large- 
scale dynamo. In order to keep matters simple, we restrict 
ourselves to the case of an isothermal equation of state with 
constant sound speed c s . Our goal is then to analyze the 
appearance of the resulting magnetic field above the surface 
of the dynamo and to study also the effects on the dynamo 
itself. 



3. The model 

The idea is to combine the evolution equations for the dy- 
namo interior with those of the region above by simply 
turning off those terms that are not to be included in the 
upper part of the domain. We do this with error function 
profiles of the form 



6 w (z) = \ (l - erf^) , 



(3) 



where w is the width of the transition. Thus, the momentum 
equation is assumed to take the form 



^ = 6 w (z) (-V/i + /) + J x B/p - 



(4) 



where F v { sc = p _1 V • (2puS) is the viscous force, Sij = 
\(JJij + Uj y i) — \5ij V • U is the traceless rate-of-strain ten- 
sor, commas denote partial differentiation, h = c 2 In p is the 
specific pseudo-enthalpy, c s = const is the isothermal sound 
speed, and / is a forcing function that drives turbulence in 
the interior. The pseudo-enthalpy term emerges from the 
fact that for an isothermal equation of state the pressure 
is given by p = Cg/?, so the pressure gradient force is given 
by p _1 Vp = CgVlnp = V/i. The continuity equation can 
either be written in terms of h 



(5) 



where we have inserted the 9 w (z) factor to terminate the 
evolution of h in the exterior, or in terms of p, 



dp 
~dt 



-V-[e w (z)pU] 



(6) 



which serves the same purpose, but also preserves total 
mass. Most of the runs presented below are carried out 
using Equation (j5j), but comparisons using Equation (j6j) 
resulted in rather similar behavior. 

The forcing function consists of random plane helical 
transversal waves with wavenumbers that lie in a band 
around an average forcing wavenumber fcf. These waves are 
maximally helical with Vx / w fcf/, so the helicity is posi- 
tive. This type of forcing was also adopted in Brandenburg 
(2001) and many other recent papers. The profile function 
w (z) in front of the forcing term restricts the occurrence of 
turbulence mostly to the dynamo region, z < 0. The forc- 
ing amplitude is chosen such that the rms velocity in this 
region, u rms , is about 4% of the sound speed. 

We adopt non-dimensional units by measuring density 
in units of the initially constant density po, velocity in units 
of c s , and length in units of fc^ 1 , where fci = 2tt/L x is the 
minimal wavenumber in the x direction with an extent L x . 
In most of the cases reported below, the vertical extent 
is L z i < z < L Z 2 with L z \ = —L x /3 and L Z 2 = 2L x /3. 
In a few cases we shall consider larger domains where the 
domain is two or four times larger in the z direction than 
in the horizontal directions. The extent of the domain in 
the y direction is L y = L x . We adopt periodic boundary 
conditions in the x and y directions. For the velocity we 
employ stress-free boundary conditions at top and bottom, 
i.e. 



U XiZ U VjZ = U z = on z = L zU L z 



(7) 
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for the magnetic field we adopt perfect conductor boundary 
conditions at the bottom, which corresponds to 

A x = Ay = A ZiZ = onz = L zl: (8) 

and vertical- field or pseudo- vacuum conditions at the top, 
i.e., 

A XiZ = Ay iZ = A z = on z = L z2 . (9) 

Note that no mass is allowed to escape at the top. Although 
this restriction does not seem to affect the results of our 
simulations significantly, it might be useful to adopt in fu- 
ture applications outflow boundary conditions instead. 

Our model is characterized by several dimensionless 
parameters. Of particular importance is the magnetic 
Reynolds number, 

Re M u rm8 /rjk{, (10) 

where fcf is the wavenumber of the energy-carrying eddies. 
The ratio of viscosity to magnetic diffusivity is the mag- 
netic Prandtl number, Ptm = v/t). In all our simulations we 
use Ptm = 1. The typical forcing wavenumber, expressed 
in units of the box wavenumber, fcf/fci, is another impor- 
tant input parameter. In our simulations this value is 10. 
For the profile functions we take a transition width w with 
k\w = 0.1 in most of the runs, but in some cases it is 0.2. 
The magnetic field is expressed in terms of the equiparti- 
tion value, B ec[ , where B 2 q = po(pu 2 }, and the average is 
taken over the turbulence zone. We measure time in non- 
dimensional units r = tu rms k^ which is the time normalized 
to the eddy turnover time of the turbulence. As initial con- 
dition we choose a hydrostatic state, U = 0, with constant 
density p = po, and the components of the magnetic vector 
potential are random white noise in space with Gaussian 
statistics and low amplitude (10 -4 below equipartition). 

In this paper we present both direct numerical simula- 
tions and mean-field calculations. In both cases we use the 
Pencil CodeQ, which is a modular high-order code (sixth 
order in space and third-order in time) for solving a large 
range of different partial differential equations. 

4. Results 

We begin by considering first hydrodynamic and hydromag- 
netic properties of the model. In the turbulence zone the 
velocity reaches quickly a statistically steady value, while 
in the exterior it takes about 1000 turnover times before a 
statistically steady state is reached. This is seen in Figure [TJ 
where we show u rms (z) at different times. In the following 
we discuss the properties of the magnetic field that is gen- 
erated by the turbulence. 

4.1. Dynamo saturation 

Dynamo action is possible when ReM reaches a critical 
value Re M lt; that is about 0.5 in our case. The situation 
is only slightly modified compared with dynamo satura- 
tion in a periodic domain. For not too large values of ReM 
the dynamo saturates relatively swiftly, but for larger val- 
ues of ReM the magnetic field strength may decline with 
increasing value of ReM (Brandenburg & Subramanian 

1 http : //pencil-code . googlecode . com 
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Fig. 1. Vertical dependence of the rms velocity at different 
times. Note the drop of u rms (z) from the turbulence zone 
to the exterior by a factor of about 3-5 The inset shows 
^rms(^) hi the exterior at different times. 

2005). An example of the saturation behavior is shown in 
Figure [2j where ReM ~ 3.4. The initial saturation phase 
(100 < r < 500) is suggestive of the resistively slow satura- 
tion found for periodic domains (Brandenburg 2001), but 
then the field declines somewhat. Such a decline is not nor- 
mally seen in periodic domains, but is typical of dynamo 
action in domains with open boundaries or an external halo 
(see Fig. 5 of Hubbard & Brandenburg 2010a). The field 
strength is about 78% of the equipartition field strength, 

In all cases the magnetic field is strongest in the tur- 
bulence zone, but it always shows a systematic variation 
in one of the two horizontal directions. It is a matter of 
chance whether this variation is in the x or in the y di- 
rection. Comparison of different runs shows that both di- 
rections are about equally possible (see below). Also, the 
magnetic field pattern shows sometimes a slow horizontal 
migration, but this too seems to be a matter of chance, as 
will be discussed below. 

4.2. Arcade formation 

After some time the magnetic field extends well into the 
exterior regions where it tends to produce an arcade-like 
structure, as seen in Figures [3] and HJ The arcade opens up 
in the middle above the line where the vertical field com- 
ponent vanishes at the surface. This leads to the formation 
of anti-aligned field lines with a current sheet in the mid- 
dle; see Figures [3] and HJ The dynamical evolution is seen 
clearly in a sequence of field line images in Figure [5] where 
anti-aligned vertical field lines reconnect above the neutral 
line and form a closed arch with plasmoid ejection above. 
This arch then changes its connectivity at the foot points 
in the sideways direction (here the y direction) , making the 
field lines bulge upward to produce a new reconnection site 
with anti-aligned field lines some distance above the sur- 
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Fig. 2. Initial exponential growth and subsequent satura- 
tion behavior of the magnetic field in the interior for forced 
turbulence with dynamo action. 




the current density is shown in pink with semi-transparent 
opacity. Note the formation of a vertical current sheet above 
the arcade. 

face. Note that this sideways motion takes the form of a 
slowly propagating wave. However, it is a matter of chance 
whether this wave propagates in the positive or negative 
coordinate direction, as will be shown below in Section l4~4l 
Field line reconnection is best seen for two-dimensional 
magnetic fields, because it is then possible to compute a 
flux function whose contours correspond to field lines in the 
corresponding plane. In the present case the magnetic field 
varies only little in the x direction, so it makes sense to vi- 
sualize the field averaged in the x direction. Since the aver- 
aging commutes with the curl operator, we can also average 
the x component of the magnetic vector potential, i.e. we 
compute (A x ) x , where the second subscript indicates aver- 
aging along the x direction. This function corresponds then 




Fig. 4. Structure of magnetic field lines in the exterior, 
together with a representation of the normal component of 
the field at the interface at z = between turbulence zone 
and exterior at r = 1601. Green represents a positive and 
blue a negative value of B z . 

to the flux function of the magnetic field in the yz plane 
and averaged along the x direction. In Figure [6] we plot 
contours of (A x ) x , which correspond to poloidal field lines 
of (B) x in the yz plane. This figure shows clearly the recur- 
rent reconnection events with subsequent plasmoid ejection. 
We also compare with a color /grey scale representation of 
the x component of the x- averaged magnetic field, (B x ) x . 
Note that in the exterior the contours of (B x ) x trace nearly 
perfectly those of (A x ) x . This suggests that the x- averaged 
magnetic field has nearly maximal magnetic helicity in the 
exterior. This is also in agreement with other indicators 
that will be considered below. 

4.3. Averaged field properties 

The magnetic field is largely confined to the turbulence 
zone where it shows a periodic, nearly sinusoidal variation 
in the y direction. Away from the turbulence zone the field 
falls off, as can be seen from the vertical slice shown in 
Figure [71 Near the top boundary, some components of the 
field become stronger again, but this is probably an artifact 
of the vertical field condition employed in this particular 
case. 

In order to describe the vertical variation of the mag- 
netic field in an effective manner, it is appropriate to 
Fourier-decompose the field in the two horizontal directions 
and to define a complex- valued mean field as 

B l r(z,t) = jj l^y.z.t)^ 1 -^), (11) 

where superscripts I and m indicate a suitable Fourier 
mode. In Figure [8] we plot absolute values of the three 
components of B 01 (z,t) as a function of z for a time rep- 
resenting the final saturated state. This figure shows quite 
clearly a relatively rapid decline of \B y \ with z, while \B X \ 
and \B Z \ level off at values that are still about 40% of that 
in the turbulence zone. This suggests that our model is 
suitable to describe the evolution of magnetic fields in the 
dynamo exterior. Earlier simulations of coronal loops and 
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Fig. 7. Slice of B z through an arbitrarily chosen cross- 
section x = const. Note the periodicity with nearly sinu- 
soidal variation in the y direction, and the more nearly 
monotonous fall-off in the z direction. 
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Fig. 9. Three-dimensional visualization of the magnetic 
field viewed from above. The vertical magnetic field compo- 
nent is color-coded (yellow pointing upward and blue point- 
ing downward). Note that the field lines form a left-handed 
spiral over the scale of the domain, as expected for turbu- 
lence with positive helicity at small scales. 

4.4. Force- free versus current free 

Already the straightforward inspection of magnetic field 
lines viewed from the top suggests that the magnetic field 
is twisted and forms a left-handed spiral; see Figure [9l This 
is indeed the orientation expected for turbulence with posi- 
tive kinetic helicity, producing a negative a effect and hence 
magnetic spirals with negative helicity at the scale of the 
domain. 

We expect the magnetic field in the dynamo exterior 
to be nearly force free, i.e., we expect ((J x -B) 2 )h to be 
small compared with (-B 2 )h(«/ 2 )h- Here, (-)h denotes an xy 
average. In order to characterize the degree to which this 
is the case, we define the quantities 



v JxB 



2 {(JxB) 
Mo - 



(B 4 



,2 ((J-B) 2 )h 
^ (B 4 )h ' 



(12) 



and note that 



Fig. 8. Vertical dependence of the moduli of the compo- 
nents of B, as defined in Equation (fTTj) . 



k 2 
kjB 



1, 



(13) 



coronal heating (Gudiksen & Nordlund 2002, 2005; Peter 
et al. 2004) demonstrate that the dynamics of such fields 
is controlled by the velocity properties at their footpoints, 
which is here the interface between the turbulence zone and 
the dynamo exterior. 



k 2 
kjB 

with k 2 B = /j j q(J 2 )ii/(B 2 )b_. In Figure [TUl we show k 2 xB 
and k 2 . B as functions of z. Given that k 2 JxB /k 2 JB has values 
below 0.1, it is evident that the field is indeed nearly force 
free in the exterior. In the turbulence zone, on the other 
hand, the Lorentz force is quite significant. 

To prove the existence of a force- free structure, it is 
convenient to calculate the angle x between J and B. We 
expect x to be or tt for an ideal force- free environment. 
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Fig. 5. Time series of arcade formation and decay. Field lines are colored by their local field strength which increases 
from pink to green. The plane shows B z increasing from red (positive) to pink (negative). The normalized time r is 
giving in each panel. 
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Fig. 10. Vertical dependence of k 2 JxB /k 2 JB and kj, B /kj B . 
Note the decline of the normalized Lorentz force from more 
than 60% in the turbulence zone to less than 10% in the 
exterior. 



In Figure [TT] we show the distribution of values of \ plotted 
over the height z. One sees that, in the exterior, the angle 
X is close to 7r, but it drops to tt/2 at the very upper part, 
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Fig. 11. Two-dimensional histogram of the distribution of 
angles, p(x), where x = arccos( J • B) is the angle between 
J and B at different heights. p(x) is normalized such that 
J P(x)dx — 1- The dashed line gives the location of the 
maximum position of the distribution. 



where the normal field condition enforces that J and B are 
at right angles to each other. 

In order to demonstrate that plasmoid ejection is a re- 
current phenomenon, it is convenient to look at the evolu- 
tion of the ratio (J -B)n/ (B 2 }n versus t and z. This is done 
in Figure fT2l for L z = 6.4 and ReM = 3.4 (Run A) and in 
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Fig. 6. Time series of the formation of a plasmoid ejection. Contours of (A x ) x are shown together with a color-scale 
representation of (B x ) x ; dark blue stands for negative and red for positive values. The contours of (A x ) x correspond to 
field lines of (B) x in the yz plane. The dotted horizontal lines show the location of the surface at z = 0. 



Table 1. Summary of runs. Run A is for a cubic domain while Runs Bl to B7 are for taller domains at different values 
of ReM- Note that the nondimensional interval length Sr as well as the ejection speeds cannot be determined accurately 
for most of the runs, but the values suggests that there is no systematic dependence on the value of ReM- The column 
L lrrC gives the / and m values of the leading mode of the mean field in Equation (fTT|) and the columns 'direction' gives 
the direction of propagation of this mean field, confirming that it is a matter of chance. 



Run 


Re 


ReM 


Pr M 


k±w 


n2 /d2 
^rms/ -°eq 


~~R 2 1 R 2 


Im 


direction 


At 


Veject / Ux\ 


A 


3.4 


3.4 


1 


0.1 


0.3-0.4 


0.11 


01 


+y 


250 


0.4488 


Bl 


6.7 


6.7 


1 


0.2 


1.0-1.2 


0.52 


01 


-y 


530 


0.5464 


B2 


6.7 


13 
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0.2 
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0.45 


10 


+x 


570 


0.628 


B3 
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0.28 


01 


±y 


800 


0.5463 
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15.0 


15 
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0.29 


10 




370 


0.5236 



Figure H3] for L z = 8tt and Re M = 6.7 (Run Bl). It turns 
out that in both cases the typical speed of plasmoid ejecta is 
about 1/2 of the rms velocity of the turbulence in the inte- 
rior region. However, the time interval 5r between plasmoid 
ejections increases from w 250 to « 570 turnover times as 
we increase the kinetic and magnetic Reynolds numbers. At 
higher magnetic Reynolds numbers, the length of the in- 
terval increases to ~ 800 turnover times. A summary of all 
runs is given in Table [TJ Here we also give the I and m values 
of the leading mode of the mean field in Equation (JTTJ) and 
indicate explicitly whether the large-scale magnetic field 
varies in the x or the y direction. Indeed, both directions 
are possible, confirming that it is a matter of chance, In 



most cases the magnetic field pattern shows a slow hori- 
zontal migration, whose direction appears to be random. 
The sign in the table indicates whether the wave migrates 
in the positive or negative coordinate direction. 

4.5. Interpretation in terms of a mean-field model 

The magnetic field found in the simulations displays a clear 
large-scale structure. One may have expected that the mag- 
netic field varies mainly in the z direction, because this is 
also the direction in which the various profile functions vary. 
However, this is not the case. Instead, the main variation 
is in one of the horizontal directions (see Figure [7J where 
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simulations. The excitation condition for the dynamo can 
be quantified in terms of a dynamo number that we define 
here as 



(16) 



where 77x0 = Vto + V is the total magnetic diffusivity and 
ki was defined in Section [3] as the smallest horizontal 
wavenumber that fits into the domain. If the turbulence 
zone were homogeneous and periodic in the z direction, the 
critical value of C a is unity, but now the domain is open in 
the z direction, so one expects the dynamo to be harder to 
excite. In the models presented below we therefore adopt 
the value C a = 2.5, which is also compatible with esti- 
mates of the critical value from the simulations, if we write 
a w u rms /3 and rj t « u rm8 /3kf. 

Next, we consider the nonlinear regime by employing 
the dynamical quenching model (Kleeorin & Ruzmaikin 
1982; Blackman & Brandenburg 2002). We assume that 
a = (%k + «m, where now ax(z) = aoO w (z) is the kinetic 
a effect profile used earlier in the solution to the kinematic 
equations, and ctM is the solution to the dynamical quench- 
ing equation, 



dt 



= -277 t fcj 



SB 

B 2 

eq 



2r]k 2 a M , 



(17) 



where £ = aB — rjt^oJ is the mean electromotive force. 
Note that we have here ignored the possibility of magnetic 
helicity fluxes that must become important at larger values 
of ReM- 

In Figure [TH we compare the resulting magnetic field 
geometry with that of the direct simulations. In both cases 
the horizontal variation of the field is similar. However, in 
the direct simulations the field extends more freely into the 
exterior. This is probably caused by a vertical upward flow 
that appears to be driven by the magnetic field. In Figure[T5l 
we see the vertical dependence of an upward flow, which 
soon reaches a statistically steady state. 



the field varies mainly in the y direction). The magnetic 
field does of course also vary in_t_he z direction, but this 
happens without sign change in B z . Above the surface at 
z = 0, the field gradually decays and retains only rather 
smooth variations. 

In order to compare with dynamo theory, we now solve 
the usual set of mean-field equations for the mean magnetic 
vector potential A, where B = V x A is the mean magnetic 
field and J = V x B/fiQ is the mean current density, 

D~A _ 

— =aB-(j) t + rj)fi J. (14) 

We recall that rj = const is the microscopic magnetic dif- 
fusivity, which is not negligible, even though it is usu- 
ally much smaller than rj t . We consider first the kinematic 
regime where a and rjt are independent of B. In order to 
account for the fact that there is no turbulence above the 
turbulence zone, we adopt the profile w (z) for a and 77 t , 
i.e., we write 

a(z) = a w (z), rjt(z) = rj to w (z), (15) 

where ao and rj t o are constants, and for w we take the same 
value as for the other profile functions used in the direct 



5. Conclusions 

Our first results are promising in that the dynamics of the 
magnetic field in the exterior is indeed found to mimic open 
boundary conditions at the interface between the turbu- 
lence zone and the exterior at z = 0. In particular, it turns 
out that a twisted magnetic field generated by a helical dy- 
namo beneath the surface is able to produce flux emergence 
in ways that are reminiscent of that found in the Sun. 

Some of the important questions that still remain open 
include the presence and magnitude of magnetic helicity 
fluxes. In the present model we expect there to be diffusive 
magnetic helicity fluxes associated with the vertical gradi- 
ent of magnetic helicity density. A related question concerns 
the dependence on the magnetic Reynolds number. One ex- 
pects magnetic helicity fluxes to become more important at 
large values of ReM • 

One of the future extensions of this model includes the 
addition of shear. In that case one might expect there to be 
strong magnetic helicity fluxes associated with the Vishniac 
& Cho (2001) mechanism that may transports magnetic 
helicity along the lines of constant shear, although more 
recent considerations now cast doubt on this possibility 
(Hubbard & Brandenburg 2010b). One should also keep in 
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Fig. 14. Comparing the average in the x direction of the 
magnetic field with the 2D mean-field model. The struc- 
tures are very similar. 

mind that the magnetic field cannot really be expected to 
be fully helical, as was assumed here in order to promote 
large-scale dynamo action under relatively simple condi- 
tions. Reducing the degree of helicity makes the dynamo 
harder to excite. On the other hand, shear helps to lower 
the excitation conditions, making it again feasible to ob- 
tain large-scale dynamo action even at low relative helicity 
of the driving. Another promising extension would be to 
move to a more global geometry, including the effects of 
curvature and gravity. This would allow for the emergence 
of a Parker-like wind that turns into a supersonic flow at 
sufficiently large radii. This would also facilitate the re- 
moval of magnetic field through the sonic point. 
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